The Effect of Static Disorder and Reactant Segregation on the A + B 
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We derive the long-time behavior of the A + B —* reaction in two dimensions, finding a uni- 
versal exponent and prefactor in the absence of disorder. Sufficiently singular disorder leads to a 
(sub) diffusion- limited reaction and a continuously variable decay exponent. Pattern matching be- 
tween the reactant segregation and the disorder is not strong enough to affect the long-time decay. 
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I. INTRODUCTION 

The A + B — > reaction shows quite different behavior 
from the A + A — > reaction because of the distinguisha- 
bility of A and B reactants. Ovchinnikov and Zeldovich, 
Toussaint and Wilczek, and Kang and Redner first no- 
ticed the diffusion-controlled kinetics that can arise from 
reactant segregation Heuristic arguments predict 

that the concentration decays from an initially random 
configuration as cpjf) ~ Cdt~ d / 4 in d dimensions for 
d < 4 This scaling arises essentially from the 

competition between diffusion and randomness in the ini- 
tially random reactant densities. Simulation studies have 
confirmed these analytical results (see, for example, ||). 
The results have been extended to the case of fractal me- 
dia (see, for example, J(|[7j ). It turns out, however, that 
two dimensions is the upper critical dimension for this 
system, and so collective motions of the reactants be- 
come important for d < 2. By upper critical dimension, 
we mean the dimension below which perturbation theory 
diverges at long wavelengths. Standard arguments ||, 
then, suggest that mean field theories such as the stan- 
dard reaction diffusion partial differential equation will 
fail in two dimensions and below. The most important 
consequence of these collective effects, and what causes 
the standard theories to fail, is renormalization of the 
effective, or apparent, reaction rate. We will derive this 
natural renormalization of the effective rate from renor- 
malization group arguments. Disorder in the medium 
can also influence the reaction kinetics. The most im- 
portant effect is renormalization of the effective diffusion 
constant. Anomalous diffusion can occur if the quenched 
defect potential is sufficiently long ranged. 

The reaction A+B — > with Poissonian random initial 
conditions has recently been analyzed via field theoretic 
techniques gj. By placing bounds on the exact, stochas- 
tic partial differential equations (PDEs) for the concen- 
tration profiles, it was shown that (c^(t)) ~ Cdt~ d / 4 
for 2 < d < 4, and the coefficient Cd was determined. 
An approximate mean field theory was shown to apply 
for long times for d > 2. Rigorous bounds on the decay 
have been derived in all dimensions for the case of infinite 



action in two dimensions in the presence of correlated 
disorder. We describe the reaction and our motivation 
for this study in section II. We express the reaction in 
terms of a master equation. In section III we map the 
master equation onto a field theory using the coherent 
state representation §,[n]|l2). In section IV we derive 
an exact, stochastic partial differential equation that can 
be used to analyze this reaction system. We state the 
most natural form of mean field theory and show that it 
applies for both short and long times. In section V we 
analyze the long-time behavior of the field theory using 
renormalization group theory, finding segregation of re- 
actants into isolated domains. We match the RG flows 
to a result that is asymptotically exact in the diffusion- 
limited, matching limit. Through this procedure we are 
able to determine the universal exponent and prefactor, 
Cd, for the concentration decay. We confirm our results 
with numerical solutions of the exact, stochastic PDEs 
for the concentration profile. We then turn to the effect 
of disorder on the segregation. For sufficiently singular 
disorder, the decay exponent changes. We calculate the 
new exponent, which comes solely from the fact that the 
reaction is (sub)diffusion limited. We conclude with a 
discussion of our results in section VI. 
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The mathematical problem we consider is one where 
two distinct types of reactants, A and B, diffuse on a two- 
dimensional square lattice. When the reactants occupy 
the same lattice site, they react with reaction rate A. The 
reactants do not interact other than through the reaction. 
The reactants do, however, move in a potential field, v, 
created by disorder. The energy of reactant A at position 
x is i>(x), and the energy of reactant B at position x is 
—v(x). We take the correlations of the potential field to 
be long-ranged. Specifically, we consider the correlation 
function whose Fourier transform is given by x vv (k) — 
J <ixexp(ik • x)x OT (x) = 7/fc 2 . It turns out than any 
disorder with a correlation function less singular than 
this at long wavelengths is technically irrelevant. Such 
disorder does not lead to any interesting, new scaling 
behavior. 
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Given this description, we can state the master equa- 
tion that governs changes in the configuration of the lat- 
tice PJl2| . The master equation relates how the probabil- 
ity, P, of a given configuration of particles on the lattice 
changes with time: 

dP({mi},{ni},t) 
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Here rrii is the number of A particles on site i, and rii 
is the number of B particles on site i. The summation 
over i is over all sites on the lattice, and the summation 
over j is over the nearest neighbors of site i. The lattice 
spacing is given by Ar. The diffusive transition matrix 
for hopping from site i to a nearest neighbor site j is given 
by TA = [1 + 0( Vi - Vj )/2] and T? = [1 - (3(v t - Vj )/2]. 
The inverse temperature is given by (3 = l/(fcsT). 

The formal motivation for our study is theoretical in- 
vestigation of the effects of disorder on the A + B — > 
reaction in two dimensions. Our model, however, could 
be considered roughly to represent a reaction between 
ions on the surface of a crystalline ionic lattice. The cu- 
bic substrate lattice has dislocation line pairs, which form 
line vacancies or line interstitials. These defects are static 
and generate a random, quenched electrostatic potential 
on the surface. Each line charge produces an electrostatic 
field on the surface of the crystal that grows logarith- 
mically with distance from the charge. Charge neutral- 
ity, however, implies that the two-dimensional density- 
density correlation function of the defects should vanish 
as k 2 for small k. Since the potential is given by the con- 
volution of the defect density with Coulomb's law, the 
appropriate form of the potential-potential correlation 
function on the surface of the crystal is that assumed 
above. Moreover, the statistics of the potential field are 
very likely Gaussian at small k |T^]. On the surface of 
the crystal, the reaction A + B — > occurs, where A and 
B are ions of opposite charge. The ions interact with a 
1/r Coulomb potential, which turns out to be irrelevant 
in two dimensions. Technically, this means that the 1/r 
interaction will not affect any flow equations produced 
by renormalization group theory. This interaction can 
only affect the "matching limit." In order to avoid any 
significant effects in the matching limit, and to facilitate 
segregation of the reactants, we may want to include a 



low density of counter ions. The density of counter ions 
should be lower than that of the defects, so as not to 
screen out the effects of the disorder as well. We will not 
consider these details further, however, and we will be 
content to analyze the mathematical problem defined by 
Eq. (Q) of reactants A and B interacting with the random 
potential but not with each other. 



III. MAPPING TO THE FIELD THEORY 

We will write the reaction in terms of a field theory 
P, ^2|Jl4| . For initial conditions, we will assume that 
the A and B particles are placed at random on the lat- 
tice. The initial concentrations for A and B at any 
given site will, therefore, be Poissonian random num- 
bers. So as to reach the most interesting scaling limit, 
we take the initial average densities to be the same, 
(ca(x, 0)) = (cb(x, 0)) = no. For simplicity, we assume 
Z?A = Db = D. No adsorption or unimolecular desorp- 
tion is allowed. The random potential is incorporated 
with the replica trick ]l5| ], using N replicas of the orig- 
inal problem. The concentrations averaged over initial 
conditions are given by 



(c A (x,i)) = hm (a(x,t)) 
(c B (x,<)) = lim (b(x,t)) , 



(2) 



where the average on the right hand side is taken with 
respect to exp(— S). In the field-theoretic formulation, 
the random initial conditions have already been averaged 
over. The fields a, a, b, b remain to be integrated over in 
the field theory. Before averaging, they have no physical 
meaning and may assume complex values. The action is 
given by S = So + Si + 5 2 + S3: 

S = j d d x j 1 dta a (x,t) [d t - DV 2 + 5(t)] a a (x,£) 



+ / d d x / dtb a (x,t) [d t - DV 2 + S(t)] b a (x,t) 



51 = J d x J dt Xia a (x,t)a a (x,t)b a (x,t) 

+A 2 5 Q (x, t)a a (x, t)b a (x, t) 
+A 3 a Q (x, t)a a (x, t)b a (x, t)b a (x, t) 

5 2 = -n J d d x [5a(x,0) +b a (x,0)] 

3 2 D 2 f 

5 3 = ^— / dhdt 2 



'kik 2 k3k4 
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x a ai (ki,t 1 )a 0l (k 2 ,ti) - 6 Q2 (ki,ii)6 Q2 (k 2 ,ii) 

x a a3 (k3,t 2 )a Q 3(k4,i 2 ) - & Q4 (k 3 , £ 2 )& Q4 (k 4 , t 2 ) 
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xki • (ki + k 2 )k 3 • (ki + k 2 )x w (|ki + k 2 |) • (3) 

Summation is implied over replica indices. The nota- 
tion J k stands for J d d k/(27r) rf . The upper time limit 
in the action is arbitrary as long as tf > t. The term 
So represents simple diffusion, without an external po- 
tential. The delta function enforces the initial condition 
on the free field propagator. The term Si comes from 
the reaction terms. The parameters are related to the 
conventional reaction rate, A, by \ = A. The term S2 
comes from averaging over the Poissonian random initial 
condition. The term S3 comes from averaging the theory 
over the random potential with the replica trick. The 
potential is assumed to be Gaussian, with zero mean and 
correlation function \w (?") ■ We take 



Xw(k) = 7/fc 2 



(4) 



This correlation function leads to anomalous diffusion in 
two dimensions. 



IV. STOCHASTIC EQUATIONS OF MOTION 

An exact expression for the concentration profiles can 
be derived by performing a Hubbard-Stratonovich trans- 
formation on (|^) and integrating out the fields a and b: 

d t a nv = DV 2 a nv + (3DV ■ (a^Vti) - \ia nv b nv 

+O/2 + im) a nv 
dtb^y = DV 2 b nv - I3DV ■ (b vv Vv) - A 2 a I)t ,6^ 
+(773 + irji)b nv 
a 7)K (x,0) = br, v (x,0) = n , (5) 

where the real, Gaussian, random field r\i has zero mean 
and variance 



(6) 



The potential field u(x) enters the equations in the ex- 
pected fashion. The physical concentration is given by 
averaging the solution over the random fields r\i . The so- 
lution before this averaging has no physical meaning — it 
is complex. Recall, also, that the random initial con- 
ditions have already been averaged over. It is for this 
reason that we have deterministic initial conditions in 
Eq. (|). 

A natural approximation for the concentration profiles 
is the standard reaction diffusion PDE with Poissonian 
random initial conditions: 



d t CAv = DV 2 c Av + j3DV ■ (ca„Vw) - Aca„cb„ 
d t CBv = DV 2 c Bv - (3DV ■ (cb„Vu) - XcxvCbv 



(7) 



Unlike the mean field theory that arises from a direct sad- 
dle point approximation to action (JsJ) || , which has uni- 
form initial conditions, these equations are valid for both 



short and long times. At intermediate times, they are 
incorrect, because they do not capture the renormaliza- 
tion of the effective rate constant. Our renormalization 
group treatment will predict how the effective reaction 
rate at length scale Are 1 scales with I. One should in- 
terpret this result as the effective, renormalized reaction 
rate observed at this length scale. With this interpreta- 
tion, an approximate means of correcting Eq. (Q) would 
be to use the running coupling A(7) in place of the origi- 
nal rate constant A. As we will see, however, the reaction 
becomes diffusion limited at long times, and so the actual 
value of A(Z) does not matter. 



FLOW EQUATIONS AND THE MATCHING 
LIMIT 



We use renormalization group theory to analyze the 
long-time behavior of the action (||) ||. The flow equa- 
tions are similar to those for the A + A — > reaction [ 14) . 
The flow equations in two dimensions, to one loop order, 
are 



din no 

dl 
din Ai 

dl 
din A 2 

dl 
din A3 

dl 

dP 2 l 
dl 
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A3 
AttD 
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= 



The dynamical exponent is given by 



z = 2 



^7 
An 



(8) 



(9) 



We note that if the sign of the attraction for the A and 
B reactants to the potential were the same, instead of 
opposite, the second term in the flow equations for A; 
would be +3/5 2 7/(4tt). 

The flow equations are integrated to a time such that 



= te- fo' z{l}dl = t 



(10) 



The matching time, to, is chosen to be relatively small 
so that mean field theory can be applied. It must be 
non-zero, however, so as to remain within the scaling 
regime. We note that for large I*, the renormalized initial 
density is large, the renormalized reaction rate is small, 
and the renormalized disorder strength remains at the 
bare value. At these short times, the diffusion is normal 
with essentially the bare diffusivity. 

We can derive the appropriate matching theory by 
considering the renormalized reaction within the original 
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master equation formulation. The Poissonian initial den- 
sities become asymptotically Gaussian since they are so 
large. We consider one lattice site and assume without 
loss of generality that ca(x, 0;Z) > cb(x, 0;/). Because 
the renormalized density is large, the reaction terms will 
dominate. We see that cb(x, t; I) will be driven to a value 
of order l/[fc(Z)cA(x, t; /)]. This will happen in a time of 
the order of l/[fc(/)cA(x, t; I)], before any significant dif- 
fusion can take place. For times t > t\, then, we will 
find that there is either a significant amount of A or B 
at each lattice site, but not both. The matching limit 
is, therefore, like the infinite reaction rate limit, as long 
as the renormalized reaction rate is not too small. The 
single-species lattice sites are the renormalized vestige 
of the long-time reactant segregation that occurs in the 
A + B — > reaction. In this infinite reaction rate limit, 
the field (p(x, t) = ca(x, t; I) — cb(x, t; I) exactly satisfies 
the diffusion equation: 



p(x,0) 



2n (l) 
Ar 2 



1/2 



AT(x) , 



(11) 



where iV(x) is a unit Gaussian variable. The potential 
has been omitted from Eq. (^IJ) because it does not affect 
the short time behavior. Furthermore, we can invert the 
definition of ip to obtain 



c A (x, t; I) = <p(x, t)H[(p(-x., t)] 
cb(x,*;Z) = — (p(x,t)H[— <p(x.,i)] 



(12) 



where H{x) is the Heavyside step function. The field ip 
is a Gaussian random variable, since it satisifies a linear 
equation and is Gaussian at t = 0. The average concen- 
tration of A and B at the matching time is, then, given 
by 



(c A {t ;l)} 



cb(*o;0) 
V(*q)) 

2vr 



1 

~ 2 

1/2 



<K*o)l) 



2no(l)G(0,2t ) 



-,1/2 



2- 



(13) 



where G(0,t) = exp(-4Dt/Ar 2 )I^(2Dt/Ar 2 )/Ar 2 is 
the two-dimensional lattice Green's function evaluated 
at x = 0. These average concentrations are indepen- 
dent of position, once we have averaged over the initial 
conditions in Eq. (JeD). For t > Ar 2 /2D, we find 



(c A (io;0) = (cb(*o;0) 



n (l) 



lTT 2 Dt Q 



1/2 



(14) 



To confirm this matching limit, we present in Figure 
1 numerical solutions of the exact stochastic PDEs (j|). 



The PDEs were integrated with an explicit Euler method 
]l6t on square 256 x 256 or 512 x 512 grids. We clearly see 
three time regimes. During the first, < t < Ar/ (ky/n ), 
the minority reactant disappears at each lattice site. 
During the second, Ar/(/c v / n ) < t < Ar 2 /D, the re- 
action is diffusion limited, but lattice effects are present. 
During the third, t > Ar 2 /D, the reaction is diffusion 
limited, and we have reached asymptotic scaling. A lim- 
ited range of the asymptotic scaling regime is shown in 
Figure 1. Simulations run to longer time would clearly 
show this scaling. Note that the range of data shown in 
Figure 1 is enough to verify Eq. ( |l4| ) in the short-time, 
matching limit, which is all we require. Furthermore, as 
shown below, the scaling of the concentration with y/n 
for all but very short times, which is clearly seen in Figure 
1, is enough to establish the exponent of the long-time 
decay. 

Alternatively, we can do direct perturbation theory on 
Eq. (g). We define ip = a — b, cf> = a + b. We use the 
reference state ciq = bo = l/(n " 1 + Ait) and perturb in 
the parameter A3. The first order corrections are 



2A 3 tdti^-t' 



„ (n^ + Aii') 2 
2 A., 



(n 1 + Ait) - ,/n 
x(n~ 1 + Ait') 



dt'G[0,2(t-t')] 



(15) 



Furthermore, ip is Gaussian to this order so we find 
(\ip\) = (2(-0 2 )/7r) 1 / 2 . To second order we find 



A, 



2(rio ^ W ^0 
x^ 1 + Ait') 2 



(16) 



In Figure 2 we plot the exact results as well as the pertur- 
bative results. We see that this perturbation theory for 
(a) is accurate until the reaction becomes diffusion lim- 
ited. In the long time regime, the perturbation theory 
fails. Interestingly, the first order result for (ip 2 ) seems 
to be exact. Asymptotically we have (ip 2 ) ~ n /(47r_Dt), 
which we have seen is correct in the long-time, diffusion- 
limited regime. But (ip 2 ) also agrees with the exact (t^ 2 ) 
for all times and for all initial densities that we have 
examined numerically. Note that the field ^(x, t) has al- 
ready been averaged over the initial conditions, and so 
it does not equal Ca(x, t) — c B (x, t) for a given instance 
of the initial conditions. This identification does hold, 
however, when each lattice site contains only A or B. 

Finally, we can also solve the standard reaction diffu- 
sion PDEs (0) with Poissonian initial conditions. Upon 
doing so, we would find that the reaction diffusion PDEs 
are accurate for short times. They are also accurate for 
long times, since they correctly capture the diffusion- 
limited nature of the reaction. Only for intermediate 



4 



times do we see a discrepancy. It is only during these 
times that the renormalization of the effective rate con- 
stant has an effect on the observed density. The discrep- 
ancy becomes smaller for larger values of no . It would not 
be noticeable for the parameters of Figure 1, for example. 

We are now in a position to combine, with confidence, 
the RG flow equations and the matching theory. We take 
to ^> Ar 2 /(2D) to be in the scaling regime. There will 
be no significant renormalization of D if to is not chosen 
too large (t = Ar 2 /D is reasonable). Using Eqs. (||), 
©, and Q we find 



(ca(*)> = (cb(*)> 



"o 



8iT 2 Dt 



1/2 ( t 

*o 



8/2 



with 



8ir 



(17) 



(18) 



VI. DISCUSSION 

We have found the exponent and prefactor for the con- 
centration decay of the A + B — > reaction at long times. 
These quantities are universal in the absence of disorder. 
Unlike the A + A — > reaction, we do not have a loga- 
rithmic correction in the upper critical dimension d c = 2 
in the absence of disorder. This is because Aj(i*) has no 
effect on the matching limit. Like the A + A — > reac- 
tion, however, the correction to the decay exponent due 
to the disorder is exactly the same correction that occurs 
to the mean-square displacement exponent This 
phenomenon is not accidental and occurs precisely be- 
cause the reaction becomes diffusion limited in the long- 
time limit. For this reason, we do not expect higher loop 
calculations of the flow equation to produce terms that 
modify /3 2 7. We have checked that the two-loop flow 
equations do not modify f3 2 j. Higher order calculations 
may modify the flow equations for A; , but this will not af- 
fect the concentration decay at long times. In fact, if the 
flow equations for Ai were unchanged by higher orders of 
/3 2 7, the matching limit would be like the infinite reac- 
tion rate limit only for f3 2 j < Air. Note that an extension 
of the RG argument gives (c A (t)) ~ ^ f) /[y/ii(8nDt) d/i ] 
for 1 < d < 4 in the absence of disorder, extending the 
results of || to one and two dimensions. 

The reactants A and B are attracted to regions of the 
potential that have opposite signs. One might expect 
that the segregation of the reactants will be biased by 
the presence of the disorder. Indeed, pattern matching 
between the reactant segregation and the potential does 
occur to some extent: when A and B have opposite at- 
traction for the potential, the effective reaction rate is de- 
creased, asymptotically reaching zero. If the attraction is 



the same, the effective reaction rate is increased, asymp- 
totically reaching a non-zero fixed point value. The dom- 
inant mechanism for the long-time decay of the reactants, 
however, is diffusion-limited decay of initial density vari- 
ations. These density variations are, in our model, not 
affected by the disorder. 
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Figure 1. Presented are numerical solutions of the exact, 
stochastic PDEs (|) for n = 10 3 , 10 4 , and 10 5 . For 
each solution we show (a) (solid), (6) (solid), (\a — b\)/2 
(short dashed), and \/((a — b) 2 ) (long dashed). We have 
used k = D = Ar = 1, /3 2 7 = in this and all figures. 

Figure 2. Presented are the perturbative results (|l5| ) 
and (|ili| ) as well as the exact results for n — 10. We 
show (a) (solid), (6) (solid), a + (a 2 ) (long dashed), 
\/{tp 2 ) (solid), and y^Vi) (short dashed). 
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